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ABSTRACT 



We report the discovery of a 3 kpc disk of few 10^ M© of dense, warm H2 in the nearby radio galaxy 3C326 N, which shows no signs 
of on-going or recent star formation and falls a factor 60 below the Schmidt-Kennicutt law. VLT/SINFONI imaging spectroscopy 
shows broad (FWHM~500 km s~^) ro-vibrational H2 lines across all of the disk, with irregular profiles and line ratios consistent 
with shocks. The ratio of turbulent and gravitational energy suggests that the gas is highly turbulent and not gravitationally bound. 
In absence of the driving by the jet, short turbulent dissipation times suggest the gas should collapse rapidly and form stars, at odds 
with the recent star-formation history. Motivated by hydrodynamic models of rapid H2 formation boosted by turbulent compression, 
we propose that the molecules formed from diffuse atomic gas in the turbulent jet cocoon. Since the gas is not self-gravitating, it 
cannot form molecular clouds or stars while the jet is active, and is likely to disperse and become atomic again after the nuclear 
activity ceases. We speculate that very low star-formation rates are to be expected under such conditions, provided that the large-scale 
turbulence sets the gas dynamics in molecular clouds. Our results illustrate that jets may create large molecular reservoirs as expected 
in 'positive feedback' scenarios of AGN-triggered star formation, but that this alone is not sufficient to trigger star formation. 
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1. Introduction 

The radio galaxy 3C326 N at z=0.1 is an excellent target to in- 
vestigate how the mechanical energy output of radio-loud AGN 
affects the surrounding gas and star formation. In spite of a 
few X 10^ Mq of dense H 2, akin to LIRGs, the star-formation 



rate (SER<0.07 M© yr"^ [Ogle et"aLl|2007| is orders of mag- 



nitude lower than ex pected from the Schmidt-Kennicutt law 
dNesvadba eraLl[2QT0. NIO hereafter). With the low SE R, an 
X-ray faint AGN (log^x = 40.6 erg s'^ IOgle et al.|2010l) , and a 
strong radio source (kinetic power log^kin > 44.6 erg s~^ , NIO) 
3C326 N provides a rare opportunity to clearly disentangle the 
effects of the jet from those of star formation and AGN radiation. 
It is one of the largest and oldest radio galaxies on the sky, with 2 
Mpc size, dynamical and spectral ages of 6-20 x 10^ yrs (Willis 
|& Strom ' 1978 ) and a radio core at 3 mm (NIO). 3C326 N is not 
an obvious member of a group or cluster, but has a similarly old 
(> 10 Gyr), massive (Msteiiar =fewxlO^^M0) companion at a 
projected distance of ~ 20 kpc. 

3C326 N is amongst the 30% of nearby 3CR radio galaxies 
which have bright line emission from warm molecular hydrogen 
observe d with Spitzer/IRS and line ratios consistent with shocks 
( [Ogle et al. 2010). The gas kinetic energy and line luminosities 
exceed the energy injection rates from star formation and AGN 
radiation, making the deposition of mechanical energy by the 
radio source the only plausible culprit (NIO). This energy injec- 
tion does not only trigger an outflow, but also heats the ambient 



~^ Based on observations carried out with the Very Large Telescope of 
ESO under program ID 385.B-0809. 
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gas through shocks, where the turbulent kinetic energy roughly 
equals the bulk kinetic energy of the outflowing gas. About half 
the molecular gas in 3C326 N is warm (T>100 K), a fraction 
lOOx greater than in star-forming galaxies, and suggesting that 
mechanical heating by the jet could be an important mode of 
suppressing the overall star formation in this galaxy. 

We present deep imaging spectroscopy of the ro-vibrational 
H2 emission lines at R=3000 and 0.7'' spatial resolution, prob- 
ing the spatially resolved gas properties and kinematics, which 
significantly enhances and complements the analysis of NIO 
which relied on unresolved Spitzer/IRS spectroscopy at R=120 
(2500 km s~^) with 5'' spatial resolution. We find a kpc-sized 
rotating, turbulent disk with broad complex line profiles, sug- 
gesting that the mechanical energy deposition of the jet is not 
restricted to a small part of the disk. We argue that finding such a 
disk is difficult to reconcile with the recent star- formation history 
of 3C326 N unless the molecular gas formed in response to the 
turbulence created by the radio jet, indicating a close symbiosis 
between jet and molecular gas, but also suggesting that boost- 
ing the formation of molecular gas in turbulent jet cocoons is 
not sufficient to trigger star formation as often assumed in mod- 
els of positive AGN feedback. Throughout the paper we adopt a 
Ho =70 km s~^ Qm =0.3, =0.7 cosmology. 



2. The molecular gas in 3C326 N 

Data were obtained with th e near-infrared imaging spectrograph 
SINEONI ( [Bonnet et aLl[2004) at the Very Large Telescope 
of ESO with a total of 15,000 seconds of on-source ob- 
serving time in the K-band at a seeing-limited resolution of 
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Fig. 1. Integrated K-band spectrum of 3C326 N. The properties 
of labeled lines are listed in Table [T] 



FWHM=0.75''x0.65''. O ur data reduction and calibration m eth- 
ods are described by, e.g Nesvadba et al. ( 2QQ8| 2011a|b| ). The 
integrated spectrum of 3C326 N (Fig.[lJ) is dominated by the ro- 
vibrational H2 lines H2 1-0 S(l) to S(4) (S(l)-S(4) hereafter). 
Paa is very weak in comparison, Paa/S(3)=0.5 (Table[T]). 

Several arguments suggest this gas is heated by shocks, sim- 
ilar to the pure-rotational H2 and optical lines (NIO). First, 
S(3)/Paa and S(l)/Br7 ratios are factors 10 greater than the typ- 
ical S(l)/Br 7^ 0.1 — 1.5 in star -forming galaxies dominated by 
UV heating ( Puxley et al.|1990 ). Since we do not detect Bry, we 
estimate S(l)/Br7=18-22 from the Paa flux and a line decrement 
Paa/Br7=12-14, which exceeds ratios in star-forming galaxies 
by a factor 10. Second, for star formation we expect ratios of the 
H2 2-1/ H2 1-0 S(3) lines factors of a few greater than the 3 a 
upper limit of 0.12 for S(l) and S(3) ( Le Petit et aL][2QQ6 



see 



Herrera et al.|201 l| for a similar analysis). A third d iagnostics is 



the [Q I]A6300/Ha ratio as a function of S(l)/Br7 ( |Mouri et aL 
19891 ). For the [OI]A6300/Ha ratio of NIO the S(l)/Br7 ratio in 



3C326 N exceeds that expected from star formation by a factor 
10. Forth, the H2 line luminosity of S(l)-S(3), 1 x 10"^^ erg s'^ 
exceeds the X-ray luminosity of the AGN and the mechanical 
luminosity of star formation (NIO), leaving shock heating by the 
jet as the only plausible culprit. 

Paa and H2 morphologies are disk-like and extend over 
2.0''x 1.1'' (3.2 kpc X 1.8 kpc) at PA^^,^/,=0±20^ All H2 lines 
and Paa are broad, FWHM--400-700 km s~^ down to the limit 
of our spatial resolution. They have irregular profiles (Fig. [2]) that 
cannot be fit with single Gaussians. To map the kinematics, we 
therefore measure the width, Wp, and central wavelength, Cp, of 
each line at p =0.2, 0.5, 0.7, and .9 x the line core extr acted 
from smaU apertures of 0.4'' x 0.4'' ( [Heckman et al.|198li . W50 



corresponds to the FWHM, and C90 approximates the line core 
For S(3) we also measure the asymmetry parameter A20, i.e., the 
shift of C20 compared to C90 (Fig. 2 of He ckman et al.||1981 



S(l) and Paa have telluric absorption features along their wings 
which makes measuring A20 difficult, but their C90 and Wso give 
results consistent with S(3). 

Velocities increase monotonically from SSW to NNE along 
PA=345°±15° (Fig. [3]) common to all percentiles and all lines. 
Lines are broadest NE of the nucleus (Wso^hi <650 km s~^ 
Wso^Pa <750 km s~^), and more narrow in the south (300-400 
km s~^ for all lines). A Fourier analysis of non-circularities in 



the velocity map with Kinemetry ( [Krajnovic et al.|[2QQ6 ) sug- 
gests the velocity gradient is indistinguishable from rotation with 
P^kin = 340°, inclination / = 35° and a deprojected circular ve- 
locity of Vc = 1/2 Avproj esc/ = 290 km s~^. This would im- 
ply a dynamical mass Mj^„ = vl R/G = 3x 10^^ Mq somewhat 
larger than the enclosed stellar mass of Mstei = 2 x 10^^ Mq. 
(The enclosed stellar mass is derived from the total stellar mass 
of 3 X 10^^ M0 (NIO) for a mass profile consistent with a Hubble 
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Fig. 2. Emission-line profile of H2 1-0 S(3) extracted from a 
0.6" X 0.6" aperture on the Northern disk (black line) fit with 3 
Gaussian components (light blue line). The fit residual is shown 
below. 



law.) The non-Gaussian line profiles suggest non-circular veloc- 
ity components which are probably smoothed out by the low spa- 
tial resolution. This could affect the velocity gradient. 

The FWHMs are larger than the velocity gradient, hence the 
broad lines cannot be artifacts of the velocity gradient and beam 
smearing, in which case most of the broadening should be near 
the nucleus, contrary to Fig. [3] The analysis of a toy data cube 
with constant FWHM=400 km s~^ and a gradient as observed 
suggests that beam smearing accounts for < 15% of the FWHM 
(similar to our Icr observational errors). We also used a toy data 
cube to quantify the contribution of rotation to the CO(l-O) line 
measured by NIO which has FWHM=350 km s~^. We find a 
minimal intrinsic line width of 250 km s~^ assuming that the 
CO has the same velocity gradient as H2 and Paa and uniform 
surface brightness. 

The small ratios ofvc/cr = 1.7 — 0.9 (with velocity disper- 
sion G = FWHM/2355 = 170 - 270 km s-^ Fig. ^ imply 
an ellipsoidal configuration rather than a thin disk. Broad lines 
across all of the disk show that the mechanical energy of the 
AGN affects the gas globally. This differs strongly from the sce- 
nario of jet-cloud interactions which are confined to gas in small 
areas along the jet axis and is more akin to recent hydrodynamic 
models of jet cocoons where radio-emitting plasma permeates 
inhomogeneous gas disks along relatively low-density channels 
( Sutherland & BickneUpOOTl [Wagner & BickneUpoTT] ). 

The H2 emission lines in 3C326 N have blue asymmetries 
which are strongest in the north (A 20 = — 150 km s~^ where tur- 
bulent velocities are also largest). H2 emission has velocities of 
up to — 1000 km s~^ from C90, well above the local escape veloc- 
ity, ^esc ^ Vc ~ 400 km s~^ at the largest H2 radii. Vesc is also 
very similar to the lowest FWHMs in the South. Fig. |2] suggests 
that of-order 10% of the emission is from gas at v > v^. Since 
the V(l-O) H2 lines only probe the warmest, but not the bulk of 
the gas, this is difficult to turn into a mass outflow rate M. If 
we assume that the pure-rotational lines, which trace most of the 
mass, have a similar profile (as recently reported for AGN with 
high-resolution Spitzer spectroscopy, |Dasyra & Combes||2011 



Guillard et al. 2011 ApJ submitted), we find 1 x 10^ M© of 
entrained dense, warm H2 in the wind. For an outflow timescale 
of 10^~^ yrs, this is consistent with or less than M estimated 
from Na D, 30 — 40 Mq yr~^ (NIO), which traces neutral (atomic 
and/or molecular) gas down to lower column densities. 
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3. Jet-disk symbiosis 

The gas in 3C326 N is in a kpc-sized, dense, thick disk. The 
broad line widths and complex profiles are not consistent with 
pure rotation, but require an additional source of turbulence 
which, following the energy arguments of NIO, can only be the 
radio source. A fraction of the line emission is from gas at veloc- 
ities above the local escape velocity. The rotational time of this 
disk. Trot = 27lR/vc = 3x10^ yrs, is shorter than the 6-20 x 10^ 
yrs age of the radio source, and comparable to the depletion time 
of the gas through the wind, T^^^^j = M/M ^ 5 x 10^ yrs. It 
is also very similar to the dissipation time of turbulent energy 
^diss=Rdisk/ ^co = 1.5 X 10^ yrs. The similar time scales raise 
two major questions. 

First, why do we find a fairly organized (albeit thick) disk, if 
the depletion time through the wind is short? The mass outflow 
rate of NIO may be uncertain to about an order of magnitude, 
but even if it were a factor 10 lower, much of the gas should 
have been removed in the 2 x 10^ yrs time span suggested by 
the oldest radio emission (Q. Finding a disk associated with 
an old radio source like 3C326 requires that energy losses are 
roughly balanced by the injection of energy, and the entrained by 
cooling gas. This would imply that a fraction of the outflowing 
gas does not escape from the galaxy, but dissipates its kinetic 
energy, stalls and rains back onto the disk in a cycle that may be 
akin to 'cold feedback' models of galaxy clusters ([Ciotti et alT] 



T99T|[Pizzoiato & Soker|2005l|Brighenti & Mathews |2006| ). 

Second, with Tdiss ^ ^rot the H2 disk cannot be stable in 
absence of the jet. The absence of a gas-rich merger suggests 
the material accumulated gradually from mass return and accre- 
tion over ^fewx 10^ yrs (NIO), longer than the age of the radio 
jet. Collisions between clouds would rapidly dissipate the tur- 
bulent kinetic energy and the angular momentum of the disk in 
about one to a few rotational times ( Verdoes Kleijn et al.|2006[ 
jPizzolato & Soker 2010) , consistent with our estimates. For the 
measured gas surface density l^gas = 250 Mq pc~^, the Schmidt- 
Kennicutt law implies SFR=4.5 Mq yr~^. Thi s is ^60 x greater 



: greater 
l.||2007| 



than the upper limit of SFR=0.07 Mq yr ^ lOg iT^ 
from Spitzer 70 fim imaging and at odds with the recent star- 
formation history of 3C326 N. To quantify recent star formation 
we repeated the population synthesis analysis of the inner 3'' 
of 3C326 N from NIO, this time adding small fractions of light 
from young (3-30 x 10^ yrs) stellar populations (YSPs) to the 
data. We find that atx SFR=4.5 Mq yr~^ stars must form for 
a few xlO^ yrs to produce detectable signatures. This is very 
short, hence, the absence of recent star formation suggests that 
the gas cannot have lost turbulent support for much longer than 





Fig. 3. (top, left to right) Maps of the line fluxes, relative veloci- 
ties, and FWHMs ofH2 l-0S(3)inkms~^ (bottom, lefttoright) 
Line asymmetry of the 20-percentile for S(l), S(3), and Paa, re- 
spectively, in km s~^. Contours mark the continuum peak. 

ally bound and be part of a turbulent multiphase medium, where 
warm H2 becomes a major gas coolant ( [Guillard et al.|2009] ). 

In this scenario the molecular gas only survives the radio 
phase if it forms self-gravitating clouds while pressurized by the 
cocoon. Otherwise, with decreasing pressure H2 formation rates 
drop, and the dynamic equilibrium between molecule formation 
and destruction is no longer maintained. The gas becomes again 
atomic. Extended reservoirs of diffuse atomic ga s are not un- 
usual in early-type galaxies (Qosterloo et al.||2007j). Dust lanes 



in rad io galaxies are often perpendicular to the jet ( |de Koff et al 



2000^, and have masses which depend on jet power (|de Ruiter 



et al.^2002) , suggesting the ISM "knows about the jet". 



4. Turbulence in 3C326 N and star formation 

To better understand how AGN and star formation in 3C326 N 
may be related, we will now argue that turbulence governs the 
gas dynamics in 3C326 N down to sizes of molecular clouds. 
In this case finding low star formation rates as observed may 
not be so surprising after all. We quantify the ratio of turbu- 
lent to gravitational energv through t he virial parameter, Oy^y = 
Sa^iinGRci I^gas)'^ (Ber toldi & Mc Kee 1992), where Z^/ and 
Rci are the mass surface density and radius of individual clouds, 
respectively, and G is the gravitational constant. Gravitationally 



a free-fall time (10^ yrs for N= 10^ cm-^) in the last 300 Myr. clouds have < 1 Hydrodynamic simulations suggest 



Unless the jet has been nearly continuously active during this 
time (longer than the oldest estimate of the jet age and unlike 
many other giant radio galaxies with signs of repeated, shorter 
jet outbreaks, |Schoenmakers et al.|2000| ), then why has this gas 
not been forming stars? 

Resolving this paradox is possible if the molecules formed 
after the onset of AGN activity, boosted by the turbulence in 
the jet cocoon. Hydrodynamic simulations suggest that turbulent 
compression can turn a diffuse inhomogeneous, largely atomic 
medium within few 10^ yrs into a warm, largely molecular 
medium ( [Glover & Mac Low||2007| ). Molecules form in rather 
dense regions (n > 1000 cm~^7"and are transported to regions 
of lower density. The total amount of H2 on macro- scales is set 
by the dynamic equilibrium of the rapid molecule formation and 
destruction on micro- scales. This gas would not be gravitation- 



that column densities created by turbulent motion cannot greatly 
exceed the mean column density (e.g., striker et_al. 2001) . We 
therefore use the measured I^gas =250 Mq pc~^ (Q- * 

The H2 line emission in 3C326 N is powered by the 
dissipation of tur bulent energy, which imp lies Lhi/^hi = 



3/2 fH2 <^ci/^d (McKee & Ostriker 2007b, where Lhi and 



Mh2 are the H2 luminosity and mass. NIO find L//2/M//2 
0.06 Lq/Mq. The correction factor fni ^ 1 is necessary because 
H2 is the dominant, but not the only gas coolant (see also Herrera 
et al. 2011, A&A accepted). We adopt ///2 = 0.5. Turbulent ve- 
locities do not scale arbitrarily with size, but approximately fol- 

as observed for molecular clouds 
1981| ), and seen in hydrody^ 



0.5 



low a power law a oc R 
in the Milky Way (e.g., [Larson 



for R=1.6 kpc. NIO found a smaller offset with a fiducial R=2.5 kpc 



namic simulations including simulations of jet cocoons ( [Krause] 
& Alexander 2007). 

Fig.|4]illustrates that for all cloud sizes < 1 kpc, the observed 
L/M implies line widths at least 10 x greater than expected for 
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Fig. 4. Turbulent velocity as a function of size. Black symbols 
show our data for CO, H2I-O lines, and Paa, respectively. The 
red dotted line is for a^ir = 1 for a constant mass-surface den- 
sity ILgas = 250 M0 pc~^ and the green hatched region shows the 
range implied by Larson's scaling law and the range of velocities 
measured on kpc scales. Black circles show the molecular clouds 
of Heyer et al. (2009 ). These points fall near the oCyir = 1 rela- 
tion for gas surface densities ^ 5 smaller than that of 3C326 N. 
The horizontal dot-dashed line shows the maximal clump size of 
-100 pc. See g]for details. 

ayir = 1 (red line in Fig. [4]), suggesting that the turbulent kinetic 
energy strongly exceeds the gravitational energy on all relevant 
scales. This includes for example the —100 pc scale at which 
disks akin to that of 3C326 N, but without additional turbulence 
from the jet, would fragment (E scala & Larson 2008 , dot-dashed 
vertical line in Fig.|4]), and also the cloud sizes of few 10s pc typ- 
ically adopted in simulations of turbulen t jet cocoons ( [Wagner & 
Bicknell|20lTt|Antonuccio-Delogu & S ilk 2008 ) 



(e.g., Krumholz & 



Fig. |4| shows also that turbulent velocities in 3C326 N are 
likely much gre ater than those in typical molecular clouds in 
the Milky Way ( [Heyer et"aL]|2009[ grey circles in Fig. |4]). This 
may hold the key to the low star-formation rates. Local density 
enhancements owing to turbulent compression may enable star 
formation in clouds with Oyfr ^ 1, but the star formation effi- 
ciency decreases rapidly with increasing 
|McKee|2005t|Padoan & Nordlund|2011| ). 

These arguments rely on the assumption that Larson's scal- 
ing relationship approximately holds for 3C326 N. Resolving 
scales of —100 pc will only become possible once ALMA 
reaches its full capabilities, however, we already see in Fig. |4] 
that Larson's scaling law (green hatched area) normalized to the 
range of velocities at kpc scales falls within the range of turbu- 
lence required to account for the observed H2 luminosity. 

5. Concluding remarks 

Our results suggest that AGN can create large reservoirs of 
dense, rapidly cooling molecular gas, as is suggested by mod- 
els of AGN-triggered star formation (e.g., Mellema et al.|2002{ 
[Fragile etal.,200 4: Silk & Norman|2009 |), but this does not per 
se seem sufficient to form stars. The turbulent energy injected 
by the jet seems capable to keep significant fractions of the gas 
warm (T>150 K) and gravitationally unbound, making star for- 
mation very inefficient. 

3C326 N is certainly an extreme example of a gas-rich radio 
galaxy with a very low star-formation efficiency, however, it is 
not unique. 30% of nearby 3CR radio galaxies have bright line 



emission from warm H2 ( Ogle et al. 2010 ), and emission-line di- 
agnostics suggest that much of this gas could be heated mechan- 
ically through shocks (NIO). This includes galaxies like 3C293, 
which have on-going star format ion, albeit not at the level ex- 
pected from their H2 reservoirs ( [Papadopoulos et aLl[2010 l. It 
will be interesting to investigate in more depth how star forma- 
tion and molecular gas kinematics are related in these cases. 
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2.2170±0.0003 


563±45 


0.7±0.05 


H2 1-0 S(3) 


1.958 


2.1339±0.0004 


594±41 


2.2±0.2 


H2 1-0 S(4) 


1.892 


2.0624±0.0004 


475±45 


0.8±0.05 


Br/ 


2.166 






<0.26 (3(7) 


H2 2-1 S(3) 


2.074 






<0.24 (3(T) 



Table 1. Line properties derived from the "stacked" spectrum, 
i.e., the integrated spectrum with the large-scale velocity gradi- 
ent removed. We also list the 3<7 upper limit for Bry, and H2 2-1 
S(3), which are not detected. 



